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We present a supersymmetric axion model in which the fermionic superpartner of 
axion, i.e. the axino, corresponds to a sterile neutrino which would accommodate 
the LSND data with atmospheric and solar neutrino oscillations. 

Current data front the atmospheric and solar neutrino experiments axe 
beautifully explained by oscillations among three active neutrino species u. 
Another data in favor of neutrino oscillation has been obtained in the LSND 
experiment □. Reconciliation of these experimental results requires three dis- 
tinct mass-squared differences, implying the existence of a sterile neutrino v s . 
the four-neutrino oscillation framework, there are two possible scenarios 
: the 2 + 2 scheme in which two pairs of close mass eigenstates are sep- 
arated by the LSND mass gap ~ 1 eV and the 3 + 1 scheme in which one 
mass is isolated from the other three by the LSND mass gap. It has been 
claimed that the LSND results can be compatible with various short-baseline 
experiments only in the context of the 2+2 scheme. However, according to 
the new LSND results □, the allowed parameter regions are shifted to smaller 
mixing n angle, thereby allowing the 3 + 1 scheme to be phenomenologically 
viable Èia. Although it can be realized in a rather limited parameter space, 
the 3+1 scheme is attractive since the fourth (sterile) neutrino can be added 
without changing the most favorable picture that the atmospheric and solar 
neutrino data are explained by the predominant — * v T and v e — + v^v T 
oscillations, respectively. In particular, the 3 + 1 scheme with the heaviest v s 
would be an interesting explanation of ali existing neutrino data. In this talk, 
we present a supersymmetric axion model with gauge-mediated supersymme- 
try (SUSY) breaking in which the fermionic superpartner of axion, i.e. the 
axino, corresponds to a sterile neutrino realizing the 3+1 scheme of four neu- 
trino oscillation u. In this model, axino can be as light as 1 eV, and a proper 
axino-neutrino mixing is induced by i?-parity violating couplings which ap- 
pear as a consequence of spontaneous U(1)pq breaking. It turns out that 
only the large angle MSW solution to the solar neutrino problem is allowed 
in this model. 

The model under consideration contains three sectors: the observable sec- 
tor, the SUSY-breaking sector, and the PQ sector. The observable sector con- 
tains the usuai quarks, leptons, and two Higgs superfields, i.e. the superficlds 
of the minimal supersymmetric standard model (MSSM) . The SUSY-breaking 
sector contains a gauge-singlet Goldstino superfield X and the gauge-charged 
messenger superfields Y, Y c as in the conventional gauge-mediated SUSY- 
breaking models Finally the PQ sector contains gauge-singlet superficlds 
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Sk (k — 1,2,3) which break U(1)pq by their vacuum expectation values 
(VEV), as well as gauge-charged superfields T,T C which have the Yukawa 
coupling with some of Sk- 

The Kàhler potential of the model can always be written as 

K = J2 + (1) 
i 

where $7 denote generic chiral superfields of the model and the ellipsis stands 
for (irrelevant) higher dimensionai operators which are suppressed by some 
powers of 1/M» where M* corresponds to the cutoff scale of the model which 
is presumed to be of order Mgut ~ IO 16 GeV. The superpotential of the 
model is given by 

W = hS 3 (SiS 2 - f PQ ) + kSìTT' + \XYY C + Wmssm + W SB (2) 

where Wmssm involves the MSSM fields, and Wsb describes SUSY breaking 
dynamics enforcing X develop a SUSY breaking VEV: (XX) = M x + 2 F X - 
This VEV generates soft masses of the MSSM fields, m so ft r Y.aFx/2nMx, as 
in the conventional gauge-mediated SUSY breaking models One can easily 
arrange the symmetries of the model, e.g. U(1)pq and an additional discrete 
symmetry, to make that Wmssm is given by 

Wmssm = y-^HiUE] + y ( g ) H 1 Q i D c j + yipH 2 QiU° + ^SfH 1 H 2 

+lLslLiH 2 + 'M-SiLiLiEi + j^SiUQ.Dl + (3) 

where the Higgs, quark and lepton superfields are in obvious notations and 
the ellipsis stands for (irrelevant) higher dimensionai operators. 

To discuss the effective action at scales below Jpq , let us define the axion 
superfield as 

A= (cj) + ia) + 6à + 6 2 F A , 

where a, <j) and à are the axion, saxion and axino, respectively. It is then 
convenient to parameterize Si and S 2 as Si — Se A 'f pQ , S 2 — Se~ A 'f pQ . As 
will be discussed later, the VEV of e 0//p « = y/Si/S 2 can be determined to be 
of order unity by SUSY-breaking effects. We then have (Si) « (S 2 ) ~ Jpq, 
and then Jpq corresponds to the axion decay Constant which would determine 
most of the low energy dynamics of axion. After integrating out the SUSY- 
breaking sector as well as the heavy fields in the PQ sector, the low energy 
effective action includes the following Kàhler potential and superpotential of 
the axion superfield A, 

KcS = f PQ {e( A+AÌ V^ + e -( A+A 'y^} + AK cS: 

W eS = (i e 2A/fpQ HiH 2 + n'ifMIfrQLiHi 

+e A ' fp * (K 3 kUL 3 El + X^UQ.Dt), (4) 
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where AK c g is A-dependent loop corrections involving the SUSY-breaking 
effects and 

= yof PQ /M*, ti = ylfpQ/Ml 
Kjk = n 3 kfpQ I M» , X' m = i ijk f PQ I M* . (5) 

The best lower bound on fpQ is from astrophysical arguments implying 
fPQ ^ IO 9 GeV 0. To accommodate the LSND data, we need the axino- 
neutrino mixing mass of order 0.1 eV. It turns out that this valuc is difficult 
to be obtained for f PQ > IO 10 GeV. We thus assume f PQ = IO 9 - IO 10 GeV 
with M* = Mgut for which /io takes an weak scale value (with appropriate 
value of yo) a and the i?-parity violating couplings Xijk, XL k are appropriately 
suppressed. It is also easy to make the coefficient \''j k of B and i?-parity 
violating operators UfDjD^ to be suppressed ^nough to avoid a too rapid 
proton decay into light gravitino and/or axino 

Low energy properties of the axion superfield crucially depends on how 
the saxion component is stabilized. One dominant contribution to the saxion 
effective potential comes from AK e s which is induced mainly by the threshold 
effects of T^Jiaving the A-dependent mass Mp = nf PQ e A/fpQ . If M T < 
M x , one fmdsE] 

167T2 Z T Z T c \ M T Ml ) ^ ' 

where Np is the number of chiral superhelds in T, Zp is the Kàhler metric 
of T, and A is a cutoff scale which is of order Mx- With Zp\ 02 g2 ~ — m^ oh , 
AK c g of (||) gives a negative-definite saxion potential 

^ (1) --^-L t l-/pQ| 2 e 2 ^-. (7) 

There is another (positive-definite) potential from the A-dependent fi- 
parameter: 

^«e^/^l/iol'd-ffiP + l^l 2 ). (8) 

With + V^ 2 \ the saxion can be stabilized at (e^^ PQ ) ~ 1 when k is of 
order 1CP 6 . Once <j) is stabilized at (e^^ PQ ) ~ 1, the resulting saxion and 
axino masses are given by 

mi « (10- 10 2 keV) 2 + Am% 

m à « (10~ 4 — 10~ 2 cV) + Am a , (9) 

where the numbers in the brackets represent the gauge-mediated contributions 
for f PQ = IO 9 - IO 10 GeV, /x Q « 300 GeV and (e^ /p «) « 1, and Am^ 
and Amg are the supergravity-mediated contributions which are of order the 
gravitino mass TO3/2 as will be discussed in the subsequent paragraph. 
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The supergravity-mediated contributions to the saxion and axino masses 
can be quite model-dependent, in particular depends on the couplings of 
light moduli in the underlying supergravity model. However they are stili 
generically of order the gravitino mass m 3 / 2 One model- independcnt 
supergravity-mediated contribution is from the auxiliary component u of the 
off-shell supergravity multiplet. In the Weyl-compensator formulation, u cor- 
responds to the _F-component of the Weyl compensator superfield: 

<5> = l + 6 2 F< 5 >, (10) 

where the scalar component of $ is normalized to be unity and the F- 
component is given by E3 

where Fj = —e~ K l 2 K IJ d{e K W^)/d&j denotes the F-component of <&/ for 
the inverse Kàhler metric K IJ which is determined by the Kàhler potential 
K of the underlying supergravity model. Note that $ is defined as a di- 
mensionless superfield, so Fq, has mass-dimension one. The above expression 
shows that F$ is generically of order m 3 / 2 . However it can be significantly 
smaller than TO3/2 in some specific models. For instance, in no-scale model 
with K = -3 In(T + T Ì > — $,$]) and dW/dT = 0, one easily finds F$ = 0. 

The Weyl-compensator contribution to the saxion and axino masses can 
be easib' read off from the super- Weyl invariant supergravity action on super- 
space E2l: 

- 3 J ée $$t e -A73 + [J d 2 e $ 3 w + h.c] (12) 

This gives the following couplings of <I> to the axion superfield: 

J éO<è&K<K = J d i 9^(A + A r ) 2 + (13) 

where K c g is the effective Kàhler potential in ([|). It is then straightforward to 
see that the Weyl compensator contributions to the saxion and axino masses 
are 

&ml = 2\F^\ 2 = 0{ml /2 ), 

Ama = F$ = 0(m 3/2 ). (14) 

In gauge-mediated SUSY breaking models i, the precise value of m 3 / 2 
depends on the details of SUSY breaking sector. However most of models give 
m 3/2 ^ 1 eV, implying that mg of Eq. (Q) is dominated by the supergravity 
contribution Ama. In this paper, we assume that Ama ~ 1 eV, so 

ma ~ 1 eV (15) 
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which would allow the axino to be a sterile neutrino for the LSND data. 
We note again that although it is generically of order 7713/2, Ama can be 
significantly smaller than m 3 / 2 when the_aupergravity Kàhler potential takes 
a particular forni, e.g. the no-scale forni lid. 

Having defined our supersymmetric axion model, it is rather straightfor- 
ward to compute the 4x4 axino-neutrino mass matrix: 



1 



(16) 



where a,8 = s, e, /z, r and v s = a with m ss = rria- The effective superpotential 
Weff in (m gives the following superpotential couplings 



d z 6 



Mo(l 



2A 
fpQ 



3A 
fpQ 



)L Z H 2 



(17) 



We will work in the field basis in which fi / i LiH2 (i = e, ji, r) in W c g are 
rotateci away by an appropriate unitary rotation of H± and Li. After this 
unitary rotation, the above superpotential couplings are changed to 



d 2 9 



2A N 

H1H2 + —, Lill2 



fpQ 



fpQ 



lcading to the axino-neutrino mass mixing 



fpQ 



0.1 



IO- 5 / V600GeV 



IO 10 GeV \ 
fpQ ) 



eV, 



(18) 



(19) 



where e, = m'ì/Mo- 



The 3x3 mass matrix of active neutrinos is induced by i?-parity violating 
couplings. At tree-level, 



(20) 



where M a denote the gaugino masses. The sneutrino VEV's (vi) are deter- 
mined by the bilinear i?-parity violations in the SUSY-breaking scalar poten- 
tial: m 2 L , H LìH\ + B' i LiH2- In our model, nonzero values of m 2 L . H and B[ at 
the weak scale arise through renormalization group evolution (RGE), mainly 
by the coupling \' i33 yb where yb is the &-quark Yukawa couplingE3. Moreovpr, 
BH1H2 arises_also through RGE which predicts a large tan/3 ~ 40 — 60 L3. 
We then find B 



m, 



10" 2 t 4 



IO" 6 



IO" 6 



eV 



(21) 



where t — ln(Mx/mj)/ln(10 3 ) for the slepton mass mr. Here we have taken 
mj ps 300 GeV and /io « 2m ; ~ which has bccn suggcstcd_to be the best param- 
eter range for correct electroweak symmetry breaking Ej. 
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Let us see how nicely ali the neutrino masses and mixing parameters are 
fitted in our framework. The analysis of Ref. [4] leads to the four parameter 
regions, R1-R4 of Table I, accommodating the LSND with short baseline 
results. In our model, Eqs. (|l5|) and ( ]Ì9| ) can easily produce the LSND mass 
eigenvalue 7714 ~ m ss = mg ~ 1 eV and also the LSND oscillation amplitude 

A — = 4£/ ^- 4 fe) 2 fe) 2 <22) 

as the four mixing elements U a 4 of the 4x4 mixing matrix U are given by 
Un f=a mi s /m ss f» 0.1 (i = e, /i,r) and U s i « 1. The masses and mixing 
of three active neutrinos can be easily analyzed by constructing the effectìve 
3x3 mass matrix given by 

< = <», - (23) 

' '"SS 

Upon ignoring the small loop corrections, this mass matrix has rank two, and 
can be written as 

mìf = mxXjXj + myYifj (24) 

where Xj and y 3 - are the unit vectors in the direction of rrii S and {vj), respec- 
tively. Remarkably, the mass scale m x « (mi S /m ss ) 2 in ss ~ IO -2 eV gives 
the right range of the atmospheric neutrino mass. Eq. (|l]) shows that m y is 
also in the range of IO -2 eV, so m cS would be able to provide the right solar 
neutrino mass unless Am 2 oi <C IO -4 eV 2 . Note from Eq. (||) that the typical 
size of e*, Xi jk ,Kjk is ar °und IO" 6 for f PQ w IO 10 GeV and M* w IO 16 GeV. 

The effective mass matrix m°f gives the mass eigenvalues 

1 



-2,3 = g 



± ^/ (m x + m y cos 2 2£) 2 + mj sin 2 2^ (25) 

and also the 3x3 mixing matrix of active neutrinos 

U 3 x3 = (z T , w T c e - x T s e , w T s g + x T c e ) ■ (26) 

where z = x x y/|x x y|, w = x x z/|x X z|, C£ = cos£ = if-^y and tan2# = 
m y sin2£/(m x + m y| e£>s2£). The Super-Kamiokande data Ii3 combined with 
the CHOOZ result Ii3 imply that C/ 2 3 « t/ 2 3 « 1/2 and J7 2 3 « 1. The 
solutions to the solar neutrino problem can have either a large mixing angle 
(LA): Ul x « C/ 2 2 « 1/2, or a small mixing angle (SA): « 1. This specify 
the first column ì T of J7|3 X 3 as 

(LA) : z « (l/x/2, -1/2, 1/2), (SA) : z « (1, 0, 0) 

up to sign ambiguities. Since x- z = 0, the pattern z ss (1, 0, 0) implies x e ss 0. 
This leads to a too small U e 4 ss m es /m ss < IO -2 , so the SA solution is not 
allowed within our model. Among various LA solutions to the solar neutrino 
problem, only the large-angle MSW solution with Am 2 oi ~ IO -4 eV 2 can be 
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naturally fitted sirice m x ss m y ~ IO -2 cV in our scheme. It is remarkable 
that fpQ « IO 10 GeV and M* ss Mqut lead to the right size of i?-parity 
violation yielding the desired values of mi S and also for the atmospheric 
and solar neutrino masses. 

To see the feasibility of our whole scheme, we scanned our parameter 
space which consists of m ss ,m,i S , X' i33 yb to reproduce the allowed LSND is- 
lands R1-R4 of Table I together with the atmospheric and solar neutrino 
parameters For RI and R4, we could find some limited parameter spaces 
which produce the corresponding oscillation parameters, however they need a 
strong alignment between x and y and also a large cancellation between m x 
and m y . On the other hand, R2 and R3 do not require any severe fine timing 
of parameters, so a sizable range of the parameter space can fit the whole 
oscillation data. 

To conclude, we have shown that the 3+1 scheme of four- neutrino os- 
cillation can be nicely obtaincd in supersymmetric axion model with gauge- 
mediated supersymmetry breaking. In this model, axino plays the role of 
sterile neutrino by having a mass ~ 1 eV and also a proper axino-neutrino 
mixing induced by i?-parity violating couplings. One interesting feature of 
the model is that only the large angle MSW solution to the solar neutrino 
problem is allowed in this model. 
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|Ami|(eV^) 


\u ei \ 




RI 
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0.11-0.13 
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Table 1. Allowcd regions for the LSND oscillation. 
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